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We propose a scheme for the calculation from the Navier- 
Stokes equations of the scaling exponents ^„ of the nth order 
correlation functions in fully developed hydrodynamic turbu- 
lence. The scheme is nonperturbative and constructed to re- 
spect the fundamental rescaling symmetry of the Euler equa- 
tion. It constitutes an infinite hierarchy of coupled equations 
that are obeyed identically with respect to scaling for any set 
of scaling exponents („■ As a consequence the scaling expo- 
nents are determined by solvability conditions and not from 
power counting. It is argued that in order to achieve such 
a formulation one must recognize that the many-point space- 
time correlation functions are not scale invariant in their time 
arguments. The assumption of full scale invariance leads un- 
avoidably to Kolmogorov exponents. It is argued that the 
determination of of all the scaling exponents requires equa- 
tions for infinitely many renormalized objects. One can how- 
ever proceed in controlled successive approximations by suc- 
cessive truncations of the infinite hierarchy of equations. Clues 
as to how to truncate without reintroducing power count- 
ing can be obtained from renormalized perturbation theory. 
To this aim we show that the fully resummed perturbation 
theory is equivalent in its contents to the exact hierarchy of 
equations obeyed by the nth order correlation functions and 
Green's function. In light of this important result we can 
safely use finite resummations to construct successive closures 
of the infinite hierarchy of equations. This paper presents the 
conceptual and technical details of the scheme. The analy- 
sis of the high-order closure procedures which do not destroy 
the rescaling symmetry and the actual calculations fWork4or 
truncated models will be presented in a forthcoming paper in 
collaboration with V. Belinicher. 



1. INTRODUCTION 

The aim of this paper is to present a general scheme 
for the calculation of the scaling exponents characterizing 
the statistical quantities that arise in the description of 
fully developed hydrodynamic turbulence. These statis- 
tical quantities are various averages computed from the 
fundamental field in hydrodynamics, the velocity field of 
the fluid. Denote the Eulerian velocity field as u(r, t) 
where r is a point in d-dimensional space (usually d — 2 
or 3) and t is the time. Statistical quantities that have 
attracted years of experimental and theoretical attention 
M-M are the structure functions of velocity differences, 



denoted as <S'„(i?) 

Sn{R) = (|u(r + R,t)-u(r,t)|") 



(1) 



where (...) stands for a suitably defined ensemble aver- 
age. It has been asserted for a long time that the struc- 
ture functions scale as a function of the separation R 
according to 



Sn{R) - R^" 



(2) 



where C„ are known as the scaling exponents of the struc- 
ture functions. It is assumed here that the separation R 
lies in the so-called "inertial range" , i.e. rj ^ R L with 
T] the inner viscous scale and L the outer integral scale of 
turbulence. One of the major questions in fundamental 
turbulence research is whether these scaling exponents 
are correctly predicted by the classical Kolmogorov 41 
theory in which C„ = n/3, or if these exponents manifest 
the phenomenon of "multiscaling" with ^„ a nonlinear 
function of n, as has been indicated by physical and nu- 
merical experiments 

In an attempt to develop a consistent theory of turbu- 
lence one may define statistical quantities that depend 
on many spatial and temporal coordinates. Defining the 
velocity difference w(r,r',i) as 



w(r,r',t) = u(r',t)-u(r,i) , 



(3) 



one considers the n-rank tensor space-time correlation 
function 

F„(ri,ri,ii;r2,r^,t2; • ■ . ;r„,r;,,t„) 
= (w(ri,ri,ti)w(r2,r^,i2) • ■ •w(r„,r;,t„)) . (4) 

The simultaneous correlation function T„ is obtained 
from F„ when ti = t2 ■ ■ ■ = tn- In statistically station- 
ary turbulence the equal time correlation function is time 
independent, and we denote it as 



T„(ri,r'i;r2,r2; . . . ;r„,r^) 
(w(ri,r'i,t)w(r2,r^,t) . . . w(r„, r^, t)) 



(5) 



One expects that when all the separations i?i = jr^ — r^| 
are in the inertial range, <C i?i <C L, the simultaneous 
correlation function is scale invariant in the sense that 



T„ (Ar 1 , Xr[ ; Ar2 , Arj ; . . . ; Ar„ , Ar^) 
= A'^"T„(ri,r'i;r2,r^; . . . ;r„,r;) , 



(6) 
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and the exponent Cn is numerically the same as the one 
appearing in Eq. (g). 

One of the major difference between the study of sta- 
tistical turbulence and other examples of anomalous scal- 
ing in physics (like critical phenomena) is that there is no 
theory for the simultaneous correlation functions (||) that 
does not involve the many time correlation functions (|^). 
Turbulence is a truly dynamical problem, and there is no 
free energy functional or a Boltzmann factor to provide a 
time-independent theory of the statistical weights. The 
theory for the simultaneous quantities involves inte- 
grals over the time variables of the many time quantities 
. One must therefore learn how to perform these time 
integrations properly. 

We propose here a point of departure from all previous 
attempts to compute the anomalous exponents from 
first principles, based on the understanding that these 
functions are not scale invariant in their time arguments; 
this allows us to build a new approach. Naively one could 
assume that the scale invariance property extends to the 
time correlation functions in the sense that 

F„ ( Ari , Ar'i , A^" 1 1 ; . . . ; Ar„ , Ar J, , A"'" i„) 
= A^"F„(ri,ri,ti; . . . ;r„,r^,i„) , (not true!) , (7) 

with a dynamical scaling exponent z„ which can be n 
dependent. We have shown that this is not the case, 
and in this paper we explain that using an adequate rep- 
resentation of the time dependence leads to a calculation 
scheme for the scaling exponents. We will reiterate that 
even in fully resummed theories, if we assume scale invari- 
ance in time we lose the availability of anomalous scaling 
in favor of K41 scaling. Full respect to the dynamical 
nature of the problem is required in order to proceed in 
understanding anomalous scaling. 

Our strategy in developing a scheme of computation 
of the scaling exponents can be described as follows: 
The first step is to transform the Navier-Stokes equa- 
tions using the Belinicher-L'vov transformation jl^, see 
Eqs.(^-(||). The resuhing velocity field is denoted W, 
and the correlation function over these velocity fields is 
denoted J^„, see Eq. (pO|). Next we construct the infinite 
hierarchy of equations that relate the rate of change of the 
nth order correlation function J^n to a space integral over 
a vertex convolved with the (n -I- l)th order correlation 
function J^n+i- This hierarchy is given by Eqs.(p6|),(pC|) 
below and it is exact. Next we explain, using the fu- 
sion rules that govern the asymptotic properties 
of the correlations functions when groups of their spatial 
arguments coalesce together, that the spatial integrals 
converge both in the UV and the IR limits. This result 
allows us to show that if Eq. (^ were true, the hierarchy 
of equations could be studied by power counting, and the 
only solution would be linear (i.e. linear in n) rather 
than anomalous scaling. We will show that Eq. {uh is not 
true, and that the many-time correlation functions must 
be "multiscaling" in their time representation. This fact 
means that the hierarchy of equations cannot be stud- 



ied by power counting. The next step is to represent 
the space and time dependent correlation functions in a 
way that exposes their multiscaling characteristics [ p^ . 
This form ( pSi) is a convenient form that makes use of the 
"multifractal" representation, which was found useful be- 
fore in the context of the simultaneous objects Q. This 
representation amounts to saddle point integrals over cor- 
relation functions that respect the exact symmetry of the 
Euler equation, i.e. the symmetry of rescaling according 
to 

r ^ Ar, w ^ A^'w, t X^'H . (8) 

Using this representation we will see that the hierarchy of 
equations for the correlation functions is satisfied in the 
sense of power counting for any value of h, cf. Eq. (|5^). 
The information about the scaling exponents is hidden in 
the coefficients of the equation but cannot be read from 
power counting. 

We are still faced with an infinite hierarchy of equa- 
tions. This is not accidental; the calculation of infinitely 
many anomalous exponents in turbulence does indeed ne- 
cessitate a calculation of infinitely many renormalized ob- 
jects. Nevertheless we want to achieve a finite, approx- 
imate calculation of lower order exponents (i.e C,n with 
for the first integer values n). To this aim we need to 
truncate the infinite hierarchy of equations. The way to 
achieve an intelligent truncation can be seen from the 
study of renormalized perturbation theory [ p^ . Pertur- 
bative theories of turbulence have fallen into some ill- 
repute since the problem of turbulence does not have any 
small parameter that guarantees convergence of pertur- 
bative expansions. To increase our confidence in this ap- 
proach, we first prove that the content of fully resummed 
perturbation theory is identical to the exact hierarchy of 
equations for n-order correlation functions and Green's 
functions. This important result, which is demonstrated 
here for the first time, allows us to proceed safely in us- 
ing partially resummed perturbation theory to guide our 
truncation of the hierarchy of exact equations. The clue 
is that we want to preserve the invariance of the the- 
ory under the rescaling (||). We outline how to do this 
in Section 8 and Appendix A of this paper. The actual 
analysis of the method of truncating while preserving the 
symmetry, and an approximate calculation of the low or- 
der exponents will be presented in a forthcoming paper 
Jisf . Here we present the conceptual steps and technical 
issues. 



2. NONPERTURBATIVE FORMULATION OF 
THE STATISTICAL THEORY OF TURBULENCE 

In order to derive anomalous exponents it is essential 
to develop a theory that is nonperturbative. In addi- 
tion, one needs to deal with the effect of the sweeping of 
small scales by large scale motions. This effect is only 
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kinematic, but it can mask the inherent time scales asso- 
ciated with nonlinear interactions, and therefore has to 
be carefully taken into account. To this aim we start this 
section with a short review of the equations of motion in 
the Belinicher-L'vov representation (subsection A), fol- 
lowed by the introduction of the statistical quantities of 
interest, which are ri-ordcr space time correlation func- 
tions and Green's functions in subsection B. In subsection 
C we derive the exact hierarchy of equations satisfie d by 
these quantities, the main results being Eqs. (^) ( |30| ) 
and iM). 



A. Equations of motion 

The analytic theory of turbulence is based on the 
Navier-Stokes equation for the Eulerian velocity field 
u(r, t). In the case of an incompressible fluid they read 

9u/ai + (u • V)u - i/V^u + Vp = f , V-u = 0, (9) 

where v is the kinematic viscosity, p is the pressure, and 
f is some forcing which maintains the flow. Since we are 
interested in incompressible flows, we project the longi- 
tudinal components out of the equations of motion. This 
is done with the help of the projection operator P which 
is formally written as P = — V^^V x Vx. Applying P 
to Eq. (^) we find 



{d/dt ~ iyV^)u + P{u ■ V)u = Pi 



(10) 



This equation has been used as a starting point for a 
field-theoretic perturbation theory, in which the nonlin- 
ear term acts as a "perturbation" on the linear part of 
the equation. Such an approach is fraught with dif- 
ficulties simply because the natural statistical objects 
which appear in the perturbation expansions 16 T9|] 
of the Navier-Stokes equation are correlation functions 
F°'P{r, r', i, t') of the velocity field u(r, t) itself: 



F^P{v,v\t,t')^{u"{v,t)uP{v\t')) 



(11) 



The problem is that the correlator F"^{r,r' ,t,t') is 
not universal, since it is dominated by contributions to 
Ua (r, t) which come from the largest scales in the fluid 
flow, and these are determined by the features of the en- 
ergy injection mechanisms. This physical fact is reflected 
in the theory as infra-red divergences that have plagued 
the development of analytic approach for decades. In- 
deed, all the early attempts to develop a consistent ana- 
lytic approach to turbulence, notably the well known the- 
ories of Wyld 1 16 1 and of Martin, Siggia and Rose pTl , 
shared this problem. One needs to transform Eq. (Mj 
such that the statistical objects that appear naturally 
will be written in terms of velocity differences which can 
be universal. We like to use the Belinicher-L'vov (BL) 
transformation for this purpose, but any other transfor- 
mation that reaches the same goal is equally acceptable. 



We will stress that the results that we obtain below do 
not depend on the details of the transformation, and we 
employ the BL transformation because it is particularly 
straightforward. The equation of motion that is obtained 
from the BL-transformation is exact, and there is no ap- 
proximation involved in this step. In terms of the Eule- 
rian velocity u(r, t) Belinicher and L'vov defined the field 
v(ro,io|i", as 



where 



v{ro,to\r,t) = u[r + p{ro,t),t] 



p(ro,i)= / dsu[ro + p(ro, s), , 

Jto 



(12) 



(13) 



Note that p(ro,i) is precisely the Lagrangian trajectory 
of a fiuid particle that is positioned at ro at time t = 
to. The field v(ro,io|r,i) is simply the Eulerian field in 
the frame of reference of a single chosen material point 
p(ro, t). The observation of Belinicher and L'vov ||ll| was 
that the variables >V(ro, to|r, r', t) defined as 

W(ro,to|r,r',t) = v(ro,io|r,i) -v(ro,to|r',t) , (14) 

exactly satisfy a Navier-Stokes-like equation in the in- 
compressible limit. Consequently one can develop a di- 
agrammatic perturbation theory in terms of these vari- 
ables [0 . The resulting theory is free of the two related 
problems that we discussed above: the (WW) correla- 
tors are universal for |r — r'| in the inertial range, and 
the theory is free of infra-red (and ultra-violet) diver- 
gences resulting from sweeping. It can be also shown 
[|ll|,0 that Kolmo gorov 1941 scaling is an order by or- 
der solution of the resulting theory. This formulation has 
two important properties, (i) The simultaneous correla- 
tors of v[ro,to|r,t] are (rg, to)-iiidependent and identical 
to the simultaneous correlators of u{r,t). The reason is 
that for stationary statistics the simultaneous correlators 
of an arbitrary number of factors of v[ro, <o|r, t] does not 
depend on t, and in particular one can take t = t^- The 
property then follows directly from Eqs. (|l^)-(|l^). (u) 
The correlators of W are closely related to the structure 
functions of u, Eq. (Q). Clearly 



5„(r-r') = (|W(ro,io|r,r',i)r) 



(15) 



As the structure functions of the Eulerian velocity differ- 
ences have for years been at the focus of experimental re- 
search, the formulation in terms of the variables W gives 
a direct link between theory and experiments. We will 
refer to the variable W(ro, io|r, r', i) as the BL velocity 
differences. We now perform the Belinicher-L'vov change 
of variables (|l^p^ together with f [r -I- p(ro, i): = 
</)[ro, to|r, t]. Using the Navier Stokes equation and the 
chain rule of differentiation we find the equation of mo- 
tion for W(ro, to|r, r', t): 



d_ 

dt 



£-KV^ + V;2) W(ro,to|r,r',0 



= P0(ro,to|r,t) -P>(ro,io|r',t). 



(16) 
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We introduced an operator L — >C(ro, to|r, r', t) as fol- 
lows: 

£(ro, io|r, r', t) = PW(ro, to|r, ro, t) • V, 

+ P'W(ro,io|r',ro,0- V; . (17) 

We remind the reader that the application of P to any 
given vector field a(r) is non local, and has the form: 



[Pa(r)]" = y dfP"^(r - v)aP{v) 



(18) 



The explicit form of the kernel can be found, for example, 
in 0] . In (H^) P and P' are projection operators which act 
on fields that depend on the corresponding coordinates r 
and r'. The equation of motion ( |l6| ) forms the basis of the 
following discussion of the statistical quantities. It has 
the tremendous advantage over the Eulerian version ( |l0| ) 
in that all the sweeping effects are removed explicitly. We 
note that the equation of motion is independent of a-nd 
from now on we drop the argument to in VV(ro, io|r, r', f). 



B. The statistical quantities 

The fundamental statistical quantities in our study 
are the many-time, many-point, "fuUy-unfused" , n-rank- 
tensor correlation function of the BL velocity differences 
>Vj = W(ro|rjr^ , tj). To simplify the notation we choose 
the following short hand notation: 



{r„ r' , t,}, X, ^ {r„ i,}, ^ >V(X,) , (19) 



^„(ro|Xi...X„) = (Wi...>V„) 



(20) 



By the term "fully unfused" we mean that all the coordi- 
nates are distinct and all the separations between them 
lie in the inertial range. In particular the 2nd order cor- 
relation function written explicitly is 



(ro|ri,ri',ti;r2,r2,i2) 
= (W"(ro|ri,r;,ii)>V^(ro|r2,r'2,i2)) 



(21) 



In addition to the n-order correlation functions the sta- 
tistical theory calls for the introduction of a similar array 
of response or Green's functions. The most familiar is the 
2nd order Green's function G°'^{vq\Xy \ x-i) defined by the 
functional derivative 



G"'^(ro|Xi;x2) 



(5W"(ro|Xi) 
<50/5(ro|x2) 



(22) 



In stationary turbulence these quantities depend on t\ — 
t-i only, and we can denote this time difference as t. 

Consider next the nonlinear Green's functions Qm.n 
which are the response of the direct product of m BL- 
velocity differences to n perturbations. In particular 



r u\x Y \ / '^[>V(ro|^i)W(ro|X2)] ' 



5i,2(ro|^i;a;2,.T;3' 



Q2.l(yf)\Xx-,X2,Xz,X^) = 



S^W{ro\Xi) 
(ro|x2)(5(?!)(ro|a;3)/ ' 
S^[W{ro\Xi)W{ro\X2)] 



(23) 



S<j){ro\x3)6<j){ra\x4) 
g3^i{ro\Xi,X2,X3;x4) 

' <5[W(ro|Xi )>V(ro|X2)W(ro|X3 )] ' 
S(f){ro\x4) 



Note that the Green's function G of Eq. ( p^ ) is Qi^i in 
this notation. 



C. Hierarchies of Equations for the Statistical 
Quantities 

1. The correlation functions 

The rate of change of the n-order correlation functions 
with respect to any of the time variables is computed as 

dti oti 



Using the equation of motion (|16|) we find 
dTn{ro\Xi, . . ■,Xn) 



Vn{ro\Xi, 



: Xn ) 



dti 

= J„(ro|Xi,...,X„) , (25) 

I?„(ro|Xi,...,X„) = {{tW)iW2 W„), 

JMX^, . . . ,X„) = i^(V? + V'i)(Wi . . . . . . W„). 

with Ci = £(ro, io|ri, r'l, t). We remember that (£yV)i 
is a nonlocal object that is quadratic in BL-velocity dif- 
ferences, cf. Eq. (p^). In writing Eq. ( p5|) we discarded 
the forcing term; it was shown before that if the forcing is 
limited to the largest scale L, its contribution is negligi- 
ble when all the separations in the correlations functions 
are much smaller than L. 

To understand the role of the various contributions in 
Eq. (25) we first note that in the limit v ^ Q the term 
Jn{'^o\Xi, . . . , Xn) vanishes. To see this note that in the 
fully unfused case this term is bounded from above by 
Ci^j7"„,p({Rj}, 0)/i?j^ijj where C is a i/-independent con- 
stant and i?min is the minimal separation between the 
coordinates. There is nothing in this quantity to balance 
V in the limit — > 0. This is of course the advantage 
of working with fully unfused quantities; we could not 
do this with the balance equation for fused correlators 
[say structure functions S'„(i?)] as the dissipative term 
approaches a finite limit when — > 0. Thus for — > 0, 
or for very large Reynolds numbers, we have 



dTn{Ya\Xi, ...,X„ 
dU 



P„(ro|Xi,...,X„) =0 . (26) 
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a. Vertices and correlation functions: 




b. Greens' functions: 




FIG. 1. The diagrammatic notation of the basic objects of 
the theory. Panel a: the vertex 7 and the correlation functions 
Tn with n = 2, 3, 4. The vertex F has not been defined yet, 
but it appears later, cf. Eq. (61). Panel b: the bare Green's 
function C/",! (thin line), and the dressed Green's functions 
Q-n,ra- Objects with only straight tails are identically zero. 



These equations, which are exact, can be written as a 
chain of equations relating the rate of change of the nth 
order correlation function to the (n + l)th order correla- 
tion function. To this aim introduce the vertex function 

7"^^(ri, r;,f) = i{[P"nri - ? ) - P"^(r; ~ f)] 

+ [P"'^ (r 1 - f ) - P^- (r'l - f )] J- ■ (27) 

With the help of this function we rewrite Eq. (|l^) in the 
form 

[^-KV? + V?,)]W"(ro|Xi) 
+ y" df7"^-(r,r',f)>V''(ro|li)W"(ro|Xi) 

= y df[P"^(ri - f) - P"^(r; - f)],^'^(ro|f,ti). (28) 

Here we used the shorthand notation 

Xi = (f,ro,ii) . (29) 

From this we can present immediately the interaction 
term I?„ in the hierarchy of equations as follows 

I?,f^-(ro|Xi,...,X„) = y df7"^'^(ri,rl,f) 

x^,':^f ■■•(ro|Xi;Xi;X2,...,X„+i) . (30) 

In order to make the structure of these equations more 
transparent we will represent them graphically using the 
diagrammatic notations introduced in Fig.^. The con- 
vention chosen here is that every velocity difference is 
represented by a wavy line of unit length, and the 2-point 
correlator as a wavy line of twice that length. The forcing 
is represented by a straight line of unit length. Thus the 
Green's function, which is the response of the velocity 
field to the forcing, is represented by a wavy line of unit 
length connected to a straight line of unit length which 
represents the forcing. This convention is generalized to 
the higher order quantities as is evident from Fig. |l|. The 
higher order quantities carry a "junction" that connects 
the appropriate number of wavy and straight lines of unit 
length each. The diagrammatic representation of the hi- 
erarchy ( p6| ) with the interaction term (^0|) is shown in 
Fig. |2| The numbers j = 1,2,3,4 represent sets of coordi- 
nates together with the corresponding vector indices 
of the correlation functions Tn. The two wavy lines des- 
ignated by "1" are connected to the empty circle 7. This 
means that both these lines carry the same set of coor- 
dinates Xi, and one is supposed to integrate over f, as is 
required by Eq. (|30|). 



2. The Green's functions 



The equation of motion for the standard Green's func- 
tion G = Qi^i of Eq. ( |2^ ) may be derived straightfor- 
wardly. Begin with Eq. ( p8| ) for >Vf and evaluate its 
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Equations for correlators: 



d ' 

-gfT Arw\i\p- + 




+ 



1 1 

2 



= 0, 



1 r^2 




= 0, 



d 





= 0.., 



1 ^ 4 



FIG. 2. The hierarchy of equation for the correlation func- 
tions. The rate of change of Tn is related to an integral over 
Tn+i according to Eqs. (29) and (33). There are two tails 
designated by 1, carrying the coordinates X\. An integration 
over f is assumed. 



functional derivative with respect to 
one gets 



After averaging 



= G(°)"'3(ri,r'i,r2)<5(ii-t2) . (31) 

The RHS of this equation displays the (zero-time) bare 
Green's function of Eq. (30), which consists of a difference 
between two transverse projection operators: 

G'(o)"/3(ri,r;,r2) = P"^{v^ -V2) - P'^^iv'^^v^). (32) 



The nonlinear term in the equation of motion (|2^) results 
in the appearance of the next-order Green's function Q2.1 
in the equation of motion (|3j). In its turn, the equation 
of motion of f/2,i will involve the next-order Green's func- 
tion Qz,i- We begin to build a hierarchy of equations for 
Qn.i similar to the hierarchy of equations for the correla- 
tion functions Tn- 

To derive this hierarchy we begin again with Eq. ( |2^ ) 
for , multiply the equation by . . . W^, and then 
evaluate the functional derivative with respect to (t)n+i- 
After averaging and discarding the viscous term the re- 
sulting equations read 



Equations for Greens' functions: 



d 

Ml 



1 1 



§(ti2) , 





= 0, 



FIG. 3. The hierarchy of equation for the Green's func- 
tions. The notation is the same as in Fig. 2. The two last 
terms on the RHS of the last equation can be obtained from 
the first one (in the same frame) by permuting 2 <-> 3 and 
2 ^ 4, as is shown. 



6 



(33) 



+ J df7"^^(ri,r;,f) 

= g(o)"'3--(i.o|Xi, . . . X„, r„+i, ti + 0),5(ii - <„+i) . 

The bare Green's function of (n, 1) order on the RHS of 
this equation are the following decomposition: 

<l"^-'^"(ro|Xi,X2...X„,r„+i,<i + 0) 
^G(°)"-(ri,r;,r„+i)^f2i-'^(ro|X2,...,X„.i) . (34) 

The diagrammatic representation of the first three equa- 
tions in the hierarchy for Qn,i are shown in FiS- The 
main difference of this hierarchy from Eqs. (^6|) is the 
presence of inhomogeneous terms on the RHS. The sec- 
ond equation in Fig. ^ is atypical in having a zero on 
the RHS. The origin of this zero is simply the fact that 
J^i = 0. 

In a similar manner we can derive additional hierar- 
chies of equations for Qn,2, Gn,3, etc. We do this starting 
from the same equation of motion (|2^) for Wi, multi- 
plying it by W2 ■ ■ ■ W„ , and then evaluating the second 
order, third order, etc. functional derivative with respect 
to the forcing. The resulting hierarchy of equations reads 



_d_ 



^^n,p(ro|^l, X2, . . . , Xn, Xn+1, 



+ y df7(ri,r'i,f) 

'>^Qn+l,p{'^Q\Xi] Xi, X2, . . . , Xn] Xn+1 ■ ■ ■ a:„+p) 
P 

= ^ G(")(ro|ri, r'l, r„+,)<5(ti - tj)e„-i,p-i(ro|X2 



(35) 



71 1 n+1 1 • 



In this equation the viscous terms have again been dis- 
carded; they are vanishingly small in the limit — > 
compared to the other terms in the equations for the 
fully unfused Green's functions. 



3. THE TEMPORAL DEPENDENCE OF THE 
STATISTICAL QUANTITIES: TEMPORAL 
MULTISCALING 

In this section we discuss the temporal dependence of 
our space-time correlation functions, culminating with a 
new representation for them, see Eq. ( p3| ) below. This 
representation will be very useful in setting up our 
scheme to compute the scaling exponents. Along the way 
we will prove the statement that if the strong assump- 
tion presented in Eq. (|^) were true, than our hierarchy 
of equations would have predicted a linear dependence of 
Cri on n. Our first step then is to show that Eq. (^ is in 
contradiction with a nonlinear dependences of C,ri on n. 



A. Locality of the Interaction Integrals 

In this subsection we discuss the locality of the inte- 
grals appearing in the hierarchic equations (|2^) and (|3^) , 
where "locality" here means that the integrals converge 
in the infrared and in the ultraviolet limits. This prop- 
erty is important, since it excludes the appearance of a 
renormalization scale as a cutoff on the spatial integrals. 
Moreover, we will show that the locality of the integrals 
in addition to the assumption of full scale invariance (7) 
imply linear scaling. Thus the only way to get anoma- 
lous scaling is to show that (7) is incorrect, as we will 
do in the next section. Fortunately, the property of lo- 
cality of the integral appearing in ( |30| ) was established 
already in previous publications. It stems from the ex- 
istence of the fusion rules which govern the asymptotic 
properties of the correlation functions when a group of 
coordinates coalesce together 13 . To control the ul- 
traviolet properties of the integrals we need to analyze 
the situation when the dummy integration coordinate f 
approaches either ro, or any of the other coordinates r;. 
The asymptotic behavior of the integrand is then deter- 
mined by the fusion rules for the fusion of two coordi- 
nates while all the rest of the coordinates remain sep- 
arated by much larger distances. The infrared limit is 
obtained when f — > cx3, and the asymptotics of the inte- 
grand are governed by fusion rules for n — 1 coordinates 
coalescing together compared to a single coordinate that 
approaches an infinite distance from all the rest. Using 
the fusion rules in this way one can prove convergence 
in a wide window of the numerical values of the scaling 
exponents around their K41 or experimental values p3[. 



B. The Consequences of Locality 

The differential equations (p^ ) can be turned readily 
into integral equations by integrating over the time vari- 
able. We argued above that the spatial integration con- 
verges in the infrared and the ultraviolet regimes. Here 
we explain the consequences of this property of locality. 
The needed material was developed in some detail in Ref. 
pO| but we repeat here all the essentials. 

In general our n-order correlation functions may de- 
pend on n different times. If Eq. (7) were correct, the 
correlation functions would be a homogeneous function 
of all the time coordinates. On the other hand, some 
of these times could be the same, and when they are all 
the same we get our simultaneous object (^ which is be- 
lieved to be homogeneous in its spatial coordinates. The 
simplest many-time case is the one in which there are 
two different times in (|o|). Chose ti = t + s for every 
i < p and ti = t for every i > p. We will denote the 
correlation function with this choice of times as T^iis), 
omitting for brevity the rest of the arguments. Introduce 
the typical decorrelation time Sn\{R) that is associated 
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with the one-time difference quantity ^„ when all 
the spatial separations are of the order of R: 

/>oo 

/ dsTi%)^siP\Ti^}iO). (36) 
Jo 

In Ref. we showed that the consequence of locality is 
that s^j^\ satisfies a p-independent scaling law: 

S„,i(i?) ~ RSn^liR)/Sn{R) oc . (37) 

There we introduced the dynamical scaling exponent z„ i 
that characterizes this time and found that 

Zn,l = 1 + Cn-1 - Cn • (38) 

Now we ask whether the same time scale also char- 
acterizes correlation functions having two or more time 
separations. Consider the three-time quantity that is ob- 
tained from J^n by choosing ti — t+si for i < p,ti — t+S2 
for p < i < p + q, and ti = t for i > p + q. We denote 
this quantity as J-^2\^1t^2)^ omitting again the rest of 

the arguments. We define the decorrelation time s^^2^ of 
this quantity by 

poo 

/ ds,ds,Tiy\s^,s,)^[si%^^]'Tiy\0,0). (39) 
Jo 

One could think naively that the decorrelation time s^^2''' 
is of the same order as (^). The calculation |0 leads 
to a different result: 

[SnAR)f - R^Sn^2{R)/SniR) CX [R'^-^f . (40) 

We see that the naive expectation implied by Eq. (^ is 
not realized. The scaling exponent of the present time is 
different from (|8|): 

Zn.2 = 1 + (C«-2 - C«)/2 . (41) 

The difference between the two scaling exponents — 
Zn,2 = Cn-1 — (Cn + Cn-2)/2. This difference is zero for 
linear scaling, meaning that in such a case the naive ex- 
pectation that the time scales are identical is correct. 
On the other hand for the situation of multiscaling the 
Hoelder inequalities require the difference to be positive. 
Accordingly, for i? ^ L we have Sn,2{R) 3> s„_i(i?). 

We can proceed to correlation functions that depend 
on m time differences. Omitting the upper indices which 
are irrelevant for the scaling exponents we denoted pC| ] 
the correlation function as ^Ti.m(si • • • Sm), and estab- 
lished the exact scaling law for its decorrelation time. 
The definition of the decorrelation time is 

/•OO 

/ ds\ . . . dSjn-F fi,m(^l • ■ ■ -^m) = [-^n^m] *^ n,m(0 • ■ • 0)- 

Jo 

(42) 



We found the dynamical scaling exponent that charac- 
terizes Sn^m when all the separations are of the order of 

Zn.m = 1 + (Cri-m - Cn)/™ , n - TO < 2. (43) 

One can see, using the Hoelder inequalities, that Zn^m 
is a nonincreasing function of m for fixed n, and in a 
multiscaling situation they are decreasing. The meaning 
is that the larger is to, the longer is the decorrelation 
time of the corresponding m+l-time correlation function, 

Sn,p{R) > Sn,q{R) for p < q . 

It is obvious now that the assumption of complete scale 
invariance (^, which is tantamount to the assertion that 
Zn^m is TO- independent, necessarily requires a linear de- 
pendence of Cn on n. An m-independent Zn^m means that 

C„-m — Cn = const X TO , linear scaling! . (44) 

The results ( ^ ) shows that in a multi-scaling situation 
our correlation functions cannot exhibit scale invariance. 

To expose the consequences in a complementary way 
that will be useful in what follows, we considered higher 
order temporal moments of the two-time correlation func- 
tions: 

/ d..'=-i^^j(.)^(;^)S^g(0). (45) 

The intuitive meaning of is^)^n \ is a k- order decorrelation 

moment of J-lf[{R, s) whose dimension is (time)*"'. The 
first order decorrelation moment is the previously defined 
decorrelation time s„ ^ • We found the scaling laws 

\^]n\ - {Sn^f - R'^Sn-km/SniR) (46) 

for k < n — 2. The procedure does not yield informa- 
tion about higher k values. We learn from the analy- 
sis of the moments that there is no single typical time 
which characterizes the s dependence of J'^l- There 
is no simple "dynamical scaling exponent" z that can 
be used to collapse the time dependence in the form 
•^r?i(*) R''" f{s/R^). Even the two-time correlation 
function is not a scale invariant object. In this respect it 
is similar to the probability distribution function of the 
velocity differences across a scale R, for which the spec- 
trum of exponents Cn is a reflection of the lack of scale 
invariance. 



C. Temporal Multiscaling Representation 

This subsection offers a convenient presentation of the 
time dependence of the correlation functions. Consider 
first the simultaneous function 7^(ri . . . rj^). Following 
the standard ideas of multifractals the simultane- 

ous function can be represented as 
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T„(ri,r;...r;) = C/" 

xT„,h(pi,pi, ■■■,p'n) ■ 



R„ 



nh+Z(h) 



(47) 



Here [/ is a typical velocity scale, and we have introduced 
the "typical scale of separation" of the set of coordinates 



n 

-E 

n ^ 



/ |2 



(48) 



Greek coordinates stand for dimensionless (rescaled) co- 
ordinates, i.e. 

-r,/i?„, =r;-/i?„, (49) 

and finally the function Z{h) is defined as 

Z{h) = 3 - V{h) . (50) 

The function Z{h) is related to the scaling exponents C,n 
via the saddle point requirement 



Cn — m\Ti[nh + Z{hy\ 



(51) 



This identification stems from the fact that the integral in 
( ^ ) is computed in the limit Rn/L via the steepest 
descent method. Neglecting logarithmic corrections one 
finds that Tn oc i?^" . 

The physical intuition behind the representation (^) is 
that there are velocity field configurations that are char- 
acterized by different scaling exponents h. For different 
orders n the main contribution comes from that value of 
h that determines the position of the saddle point in the 
integral (|^). 

There is no reason why not to extend this intuition to 
the time domain. The particular velocity configurations 
that are characterized by an exponent h also display a 
typical time scale t^^h which is written as 



(52) 



Accordingly we propose a new temporal multiscaling rep- 
resentation for the time dependent function 



^n{'^o\Xi, . . . , X„ 



dfi(h) 



Rn 



nh+Z(h) 



X-?^n,/i(ro|^l, ^2, ■ ■ • 



(53) 



The function J-n,h depends on dimensionless (rescaled) 
coordinates 



(54) 



Of 

course, we require that the fimction J-"„^/i(ro|Si, . . . S„) 



is identical to Tn^hipi, p'l, ■ • • , p'n) when its rescaled time 
arguments are all the same. 

The point of this presentation is that it reproduces 
all the scaling laws that are involved in time integrations 
and differentiations. For example, specializing to the case 
of one time difference as discussed in Section 4, we can 
perform the integral appearing in Eq. (|4^), and derive 
immediately the result in (E^). 



4. SETTING UP A CALCULATION OF THE 
SCALING EXPONENTS 

This section contains the main result of this paper, i.e. 
the equations that we propose as the starting point for 
the calculation of the scaling exponents. 

One of the main insights gained by understanding the 
temporal properties of the correlation functions is that it 
is not true that the hierarchy of equations for the correla- 
tion functions dictate classical scaling by power counting. 
The way to see this is to substitute the temporal multi- 
scaling representation ( |5^ ) in the hierarchy of equations 
(pO|), and see that power counting gives no information: 



dp{h) 



p ^ {7l+l)h + Z(h) ^ 



Rn+1 
Rn 



{n+l)h+Z{h) 







(55) 



In this equation Rn+i is the mean scale of separation 
of the 2{n + 1) space coordinates of the function Tn+i 
similarly to ( pSf ) . This set of equations can be considered 
as a linear set of functional equations for the structure 
functions in rescaled coordinates. Since this equation has 
to be valid for any value of L, and since we integrate over 
a positive measure, the equation is satisfied only if the 
terms in curly parentheses vanish. In other words, we 
will seek solutions to the equations 



_d_ 



TnM-l:---i-n)+ / dpj{pi,p2,p) 



(56) 



Rn+l 
Rn 



(n+l)h+Z(h) 



We note the important fact that this equation is invariant 
to the rescaling 



Ari 



(57) 



for any value of A and h. Accordingly, power count- 
ing leads to no result. As a consequence the informa- 
tion about the scaling exponents is obtainable only 
from the solvability conditions of this equation. In other 
words, the information is buried in coefficients rather 
than in power counting. The spatial derivative in the 
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Dyson- Wyld equations: 




1 2 1 a h 2 
^lAAAAA = AAA {^^) A/V^ 



FIG. 4. The coupled Dyson (first line) and Wyld (second 
line) equations for the Green's functions and second order 
correlator respectively. The mass operators are displayed in 
Fig. 5 



vertex on the RHS brings down the unknown function 
Z{h), and its calculation will be an integral part of the 
computation of the exponents. 

We recognize nevertheless that this set of equations 
forms an infinite hierarchy. To proceed we need to find 
intelligent ways to truncate the hierarchy without rein- 
troducing power counting. To this aim we are going to 
use renormalizcd perturbation theory, which is one of the 
best schemes available to express higher order quantities 
in terms of lower order ones. We review the needed ma- 
terial in the next section. 



5. RESUMMATIONS IN SUCCESSIVE ORDERS 



Mass operators L and o : 



+ 




B ^ 2 



FIG. 5. Exact representations of the mass operators E and 
"I> in terms of the bare vertex F, the Green's function and 
second order correlator, and the dressed three-legged vertices 
A,B,C. 



The diagrammatic representation of the Dyson- Wyld 
equations for the Green's function Qi i and the second 
order correlation function J-2 is identical to its previous 
counterpart |0, and shown in Fig. |4[ In symbols the 
Dyson equation reads now 



In this section we discuss successive resummations of 
renormalizcd perturbation theory. In doing so we will 
demonstrate that the theory generates infinitely many 
renormalizcd objects whose scaling properties are non- 
trivial. The first step is the standard line-resummation 
that produces renormalizcd two-point functions. Further 
resummations produce three-, four-, and higher point 
renormalizcd objects. We prove that this procedure gen- 
erates equations that in their fully resummed form are 
identical in content to the exact hierarchies of equations 
derived above directly from the fluid equations of mo- 
tion. Next we explain how partially resummed versions 
can be used to offer controlled approximations to the full 
calculation. 

Equation (^) can be used as a starting point for the 
development of a line-renormalized perturbation theory 
for the statistical quantities. The reader should note that 
this equation is slightly different than the one used pre- 
viously [0 for the same purpose. The difference is that 
previously one used an equation for a quantity that in 
the present notation reads W(ro|r, ro, i). Accordingly, 
one cannot read blindly the results of the previous anal- 
ysis. Nevertheless the differences are not serious, and the 
spirit is the same as before. 



+ JdraJ^' dtag[T\ro\Xi;Xa) (58) 

X J dn J^^ dtb^^-'{ro\ra,rb,ta-h)gJ^,{ro\Xb,X2) , 

where Xi, — {rh,ro,tb}. One sees that the number 1 
designates a set of coordinates Xi and a corresponding 
vector index a, the number 2 corresponds to X2 and a 
corresponding vector index (3. The meaning of "a" is 
evident from the comparison of Eq. (^8|) with the first 
line of Fig. |^. The Wyld equation for the 2nd order 
correlation function has the form 

//•OO /"OO 
dradrh / dta / dtb 
Jti Jt2 

xg'^'^{ro\Xi;Xa)g''^iro\X2;xb)^'^{ro\xa,Xb) . (59) 

The meaning of 1,2, a and b is again evident after com- 
paring the second line of Fig. ^ with Eq. ( ^9| ) . In equa- 
tion ( p8| ) the "mass operator" E is related to the "eddy 
viscosity" whereas in Eq. (^9[) the "mass operator" <& is 
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the renormalized "nonlinear" noise which arises due to 
turbulent excitations. Both these quantities are given 
as infinite series in terms of the Green's function and 
the correlator, and thus the equations are coupled. The 
diagrammatic notation for the vertex T is presented in 
Fig. [|a. Analytically Ta^j (r) is a local differential oper- 
ator of the Euler equation 



(r) = -((5qt,V/3 + (5a/3V^) 



(60) 



which is different from the non-local bare vertex 7, 
Eq. [2^. This vertex is related to T via the bare Green's 
function (ph where G^o^"^ is given by Eq. (||): 
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a/37 



(ri,r;,r2) = -GW"*(ri, r'^, r2)r^^^(r2) 



(61) 



where only the repeated tensor index S is summed upon. 
The bare Green's function in the BL-representation sat- 
isfies the equation 



d 



{Vl + Vl)]g[T'iro\X,,X2) (62) 



The series for S and 4> can be resummed exactly, as is 
shown in Fig. |^. There is a price to pay: there are three 
new objects that appear as a result of this resummation, 
known as the 'triple" vertices A, B and C. They differ 
in the number of wavy tails that connect them to the 
propagators. Note that there is no triple vertex with 
three wavy tails; such a vertex vanishes due to causality, 
as is discussed in Appendix A. Each of these vertices can 
be represented as an infinite series in terms of the same 
objects, i.e. the vertices A,B,C, and the correlator and 
propagator. The first diagrams in the infinite series for 
the triple vertices are shown in Fig. ^jb. 

In fact one is not limited to an infinite series represen- 
tation. One can also represent the triple vertices through 
exact and complete resummations of the infinite series, 
but this is at the expense of introducing yet another set 
of objects, this time of quartic nature. There are four 
different quartic vertices, with one, two, three or four 
straight "legs" . We denote them as Dm.n with m stand- 
ing for the number of wavy tails, and n for the number 
of straight tails, m + n = 4. It is natural to discuss 
the various contributions appearing in these objects as 
"skeleton" and "irreducible" and this discussion is pre- 
sented in Appendix A. 

This process goes on. We can offer an exact, fully 
resummed equation for the quartic vertices, but at the 
price of introducing quintic objects, etc, as is explained 
in detail in Appendix A. The renormalized vertices serve 
also in providing exact representations for the many- 
point different-time correlation functions (|^). As exam- 
ples we present in Figs.^],^ the exact and fully renormal- 
ized three-point correlation function J-'^ and the nonlinear 
Green's function Q2.1 in terms of the double correlator, 
the Green's function and the triple vertices. 



a. 3'rd order vertices: 

= ( ■< 1 ; =0 

ABC 



b. Series for 3'rd order vertices: 




+... 



c. Exact representation: 




+ 




FIG. 6. The 3rd order vertices. Panel a: The three different 
types of 3rd order vertices. A, B, and C. A double lined tail 
stands for either wavy or straight tail. There is no vertex with 
three wavy tails in the theory. Panel b: The diagrammatic 
representation of a dressed 3rd order vertex (A, B or C) as 
an infinite series in terms of dressed 3rd order vertices and 
propagators. There is only one type of bare vertex appearing 
in the series for the dressed vertex A. Panel c: The exact 
representation of 3rd order vertices in terms of the partly 
reducible 4th order vertex (gray square) which is defined in 
Fig. 7. 
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a. 4'th order vertices: 

X =(X .X ) . 

^4,0 = X =0 . 

b. 4'th-order vertex decomposition : 

1^4 1/^4 lA5^3 1/^4 

c. Skeleton contributions: 
d. Skeleton contributions: 

FIG. 7. The 4th order vertices. Panel a: the different types 
of 4th order vertices. A double lined tail stands for either 
wavy or straight tail, see text. Panel b, first line: decomposi- 
tion of 4th order vertices in terms of partly reducible contri- 
butions (gray square). Panel b, second line: further decom- 
position of the vertex into two-eddy reducible parts and fully 
irreducible 4the order vertices (emplty square). The latter are 
shown in Fig. 13. Panels c and d: Skeleton contributions to 
the 4th order vertices. There axe three terms in Panel d with 
three types of propagators. There is just one term in Panel 
c, the two other type of propagators require a tripple vertex 
with three wavy tails which is zero. 



3'rd order correlator: 




FIG. 8. Exact representation of the third-order correlation 
function. 



Green's functions: 




FIG. 9. Exact representation of the nonlinear Green's func- 
tions Q2,i and Gi,2- 
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4-th order correlator: 




FIG. 10. Exact representation of the 4th-order correlation 
function in terms of 4th order vertices (black squares) 
and propagators. The meaning of the 4th order vertex is 
discussed in Appendix A. The notation in this figure is that 
every contributions that is shown explicitly in a frame adds 
up to the indicated permutations enclosed in frames of the 
same line type. 



Similarly we can provide an exact representation of 
the four point correlation function .7-4 in terms of the 
same two-point correlator and Green's function and the 
quartic vertices shown in Fig. |l^. And so on; F„ can be 
represented as 2" — 1 contributions in terms of n-legged 
vertices. 

We see that we have a theory that allows exact repre- 
sentations of as many renormalized objects as we desire. 
It appears clear and elegant, but in fact is almost content- 
less. In order to endow it with predictive power we need 
to learn how to perform two tasks. Firstly, the space-time 
integrals are implicit in every diagram, see for example 
Eqs. ( p8| , p9[ ). We need to learn how to compute them. 
Secondly, we cannot deal with an infinite hierarchy, and 
we need to learn how to close the hierarchy efficiently. 
In the next section we review briefly the past history of 
attempts to solve these equations, and why they failed in 
flnding anomalous scaling. This will serve as a guidance 
for the new ideas that pave the route that is taken in the 
rest of this paper. 

6. WHY IS THE PROBLEM STILL OPEN? 

The flrst serious attempt to derive the scaling expo- 
nents of the statistical objects was due to Kraichnan 
1^,^, who introduced the "Direct-Interaction Approx- 
imation" (DIA) . In the context of our presentation this 
approximation means truncating the set of exact equa- 
tions by replacing the dressed vertex A by its bare coun- 
terpart, and by equating the dressed vertices B and C to 
their bare value which is zero. This approximation leads 
to closed equations in terms of the second order corre- 
lation functions and Green's function. Kraichnan was 
the first to understand that after an appropriate trans- 
formation of the velocity field (he chose Lagrangian, but 
Bclinicher and L'vov showed that their transformation 
amounted to the same result) the space integrals in the 
DIA converge at both ends. As a result of this "local- 
ity" of the integrals, neither an inner nor an outer scale 
appears in the bulk of the inertial interval. In addition, 
Kraichnan assumed that the space-time correlation func- 
tions are scale invariant in the sense of Eq. (^. This 
assumption, together with the property of locality, lead 
to the existence of scaling relations connecting the expo- 
nents C2 and z: 

C2 + 2z = 2 , (63) 
2C2 + ^ - 2 . (64) 

Of course, the solution is K41, i.e. C,i = z — 2/3. It 
should be understood at this point that this disappoint- 
ing result will continue to hold not just at the DIA level, 
but at any finite order in the dressing of the vertex via 
the perturbative scheme. Since it was shown that ev- 
ery diagram exhibits locality, the scaling relations are re- 
produced by any finite truncation of the line resummed 
perturbation theory. Whenever we have a set of linear 
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scaling relations, K41 can be the only solution. It is al- 
ways a solution by dimensional reasoning, and therefore 
it must solve any set of scaling relations. The linearity 
of the scaling relations is then the end of the story. 

One way out would be to assert that a full dressing of 
the vertex may result in a new scaling exponent for that 
object. We still have a linear set of two scaling relations, 
but now three exponents to determine. The extra free- 
dom allows deviations from K41. However, this is not 
the case as long as we believe that the scale invariance as 
assumed in the form Eg . (M ) indeed holds. It was shown 
by L'vov and Lebedev |23f that under that assumption 
the exponent of the vertex is unchanged from its bare 
value of —1. Accordingly, K41 seems a deeper trap than 
ever. Is there a way out? 

Our contention is that the fundamental fallacy is the 
assumption Eq. (|^) . One could think that there is a fun- 
damental problem with renormalized perturbation the- 
ory since the problem lacks a small parameter. To show 
that this is not the case we develop the theory further, 
and show that the fully renormalized equations obtained 
from resumming the diagrammatic expansions are exact 
and equivalent to hierarchical equations that stem from 
the Navier-Stokes equations. We will then argue that 
finite resummations of the perturbation theory up to a 
given order will give successive approximations which are 
useful for the solution of the fundamental equations (56). 



7. EQUIVALENCE OF THE FULLY 
RENORMALIZED DIAGRAMMATICS AND THE 
EXACT HIERARCHICAL EQUATIONS 

The student of fluid mechanics who is not familiar with 
field theoretic methods may find the diagrammatic ex- 
pansion discussed in section 5 somewhat foreboding. The 
student of field theories, on the other hand, may find 
them suspicious since the problem does not have a small 
parameter. It is important therefore for general under- 
standing, and also crucial for the further development of 
our theory, to show that the fully resummed equations 
discussed above are equivalent to hierarchies of equations 
that can be derived exactly from the Navier-Stokes equa- 
tions. 

Consider first the Dyson equations (pq). Applying to 



this equation the operator 



Kv^ + vg) 



from the 



left, and using Eq. ( |62| ) we can write, in the limit v —>■ 0^ 

J-0r?(ro|Xi;x2) =G(°)"'3(r^^r;,r2)(5(ti -t2) (65) 
oil 



dn I dtbt''''{ro\Xi;xb)gf^iro\Xt;x2) 

t2 



where 



E"^(ro|Xi;a;b) = J draG^°^'''irur[,ra) (66) 



What follows from the Dyson equation: 



5 (ti2) 




FIG. 11. The Dyson equation written in terms of the fully 
resummed vertices A and B. 



Next one might compare Fig. ITT] with the exact repre- 
sentation of the Green's function Q2.1 which is presented 
in Fig. ||. The last term in Fig. ^ is obviously a vertex 
integrated against the last contribution to C/2,1 with the 
B vertex, with a factor of 1/2. The second term is just 
the vertex integrated against the other two contributions 
with vertex A, and in total we retrieve exactly the first 
equation in the hierarchy presented in Fig. ^. 

Figure |l^ is a fiow chart of the derivation of the first 
of the hierarchy of equations for the correlation func- 
tions. The first equation, denoted by 1 over the equal 
sign, is obvious. The second one, denoted by 2, is ob- 
tained using the Wyld equation (which is the second of 
Fig. Note that the time derivative is applied only 
to the Green's function on the left. Equality 3 is a sub- 
stitution of the derivative of the Green's function from 
Fig. O. The bare Green's function (the first term in 
Fig. ^ changes the dark circle to minus an empty one, 
in accordance with Eq(|6l|)). One recognizes next that the 
pieces in the two frames designated ci are the 2nd order 
correlation function, shown in the frame C2. The first 
term shown in frame bi becomes the three contributions 
shown in frame 62 in panel b by using the exact repre- 
sentation of $ shown in Fig. ^ and relationship between 
vertices F and 7. Lastly, the equation in panel b contains 
all the contributions to the RHS, which are nothing but 
the vertex integrated against the 3rd order correlation 
function, as can be checked by comparing with Fig. |^. 
We have thus derived the first of the exact hierarchy 
of equations for the correlation functions shown as the 
first line of Fig. |^. One can proceed similarly to derive 
the second of the hierarchic equations for the correlation 
functions. The starting point would be the exact repre- 
sentation of J-3 in Fig. ^. Taking the time derivative with 
respect to ti one gets seven contributions, four of which 
involve derivative of the Green's function and three 
a derivative of ^2- For every such derivative we need to 
substitute either the three terms appearing in Fig. |l^, 
or the five terms appearing in Fig. Special attention 
should be paid to the contribution of the first term (the 
bare Green's function) in Fig. O. When it is attached 
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to one of the dressed vertices A, B or C one needs to use 
their exact representation as shown in Figs.^. In all the 
other contributions we can leave the dressed vertices A, 
B and C unchanged. At this point we can collect all the 
terms, and find 60 terms (counting each one with a fac- 
tor 1/2 once and with a factor 1 twice). One may then 
substitute the exact representation for the dark 4th order 
vertices as a sum of four contributions, like those shown 
in Fig. [T3[ b in Appendix A. This will reduce the number 
of contributions to 15, all with dark 4th order vertices. 
This is precisely the right number of contributions needed 
to reproduce the exact result which is the empty circle 
vertex 7 integrated against the 15 contributions in JF4 
shown in Fig. |l^. This completes the derivation of the 
second equation in the hierarchy shown in Fig. ^. Exactly 
the same procedure (but with fewer terms) provides the 
derivation of the second equation in the hierarchy of the 
Green's function Qn^i shown in Fig. || etc. 

This derivation can be repeated at any order, showing 
the exact correspondence of the fully resummed diagram- 
matics and the hierarchic equations. Since the proof calls 
for the introduction of additional graphical notations we 
present it in the appendix for the consideration of the 
dedicated reader. 

The summary of this section is as follows: The classical 
Dyson- Wyld equations are equivalent to the first equa- 
tions in the hierarchy of equations for the correlation and 
Green's functions. At this step of resummation the 3rd 
order objects are given only in terms of an infinite se- 
ries in the 2nd order objects. If we replace at this point 
the triple vertex by its bare value we will generate the 
Direct Interaction Approximation. In the next step of 
resummation we have 2nd and 3rd order objects in fully 
resummed form. In particular the 3rd order vertices are 
not given in terms of infinite series but they are exactly 
represented in a fully resummed form in terms of 4th or- 
der objects as shown in Fig. At this stage we can show 
agreement with the first two equations of the hierarchies. 
The 4th order objects are still presented in terms of an 
infinite series in terms 2nd and 3rd order objects, and 
one can discuss various ways of approximating the 4th 
order objects. But one can make instead the next step 
of resummation, that will yield fully resummed 4th order 
objects in terms of 5th order ones, see Fig. |l^. The 5th 
order objects are represented at this stage as an infinite 
series in terms of lower order objects. At this step we can 
recover the first three equations of the hierarchies. This 
procedure can be continued to any desired order, and at 
every n-th step of the procedure in which we have fully 
resummed n-order objects, we can recover the first n — 1 
equations of the hierarchies. It is reasonable to assume 
that by deferring the closure approximation to higher and 
higher order steps of the resummation we may find bet- 
ter and better answers for the lower order objects. In 
the next Section we will discuss several possibilities of 
making intelligent closures. 



8. CLOSURE SCHEMES 

In this section we explain why partly resummed per- 
turbation expansions are useful in implementing intelli- 
gent closures of the infinite hierarchy (p6|). We take as 
an example a closure at the level of the 3-point vertices 
A, B and C, but the comments are relevant for closures 
at higher order as well. The complete discussion belongs 
to the forthcoming paper |jl^, and here we just briefly 
sketch the important ideas, to underline why we spent so 
much effort on the development of the resummed pertur- 
bation theory in Sects. 5-7. 

The main point is that we want to select the closure 
approximation such as not to reintroduce power count- 
ing into the hierarchy ( |56| ) . Consider then the three point 
vertices A, B and C. We have for each of these vertices 
an infinite series in terms of themselves and the lower or- 
der 2-point correlation functions and Green's functions, 
see Figs.|b. Of course, taking all the infinite contribu- 
tions is impossible, and we need a criterion to choose par- 
tial series of contributions as the approximation for the 
3rd-order vertices. One traditional approximation that 
one could use is the so-called "triangular approximation" 
Ip^ ; this means that we take all the triangular diagrams 
like the first four contributions after the bare vertex in 
Fig. ||b. Each of these triangular diagrams contains three 
dressed 3-point vertices of type A, B or C, and three prop- 
agators (2-point correlator or Green's function). These 
diagrams result from the infinite (partial) resummation 
of all the "planar" diagrams which do not contain crossed 
lines of propagators. Note that the four shown triangular 
diagrams are not all the diagrams of this type, but their 
number is small (for example six in the case of vertex A). 

In fact such an approximation is not suitable for our 
needs. In the forthcoming paper [p^ we will show that 
this approximation reintroduces power counting to the 
hierarchy (|5^). On the other hand we will show that the 
approximation that is obtained by taking only the skele- 
ton contributions and neglecting the fully irreducible con- 
tributions is appropriate, since it does not reintroduce 
power counting. In other words we can use exact rep- 
resentations for the n-order fully irreducible vertices in 
terms of partly reducible n -\- 1-order vertex (as shown in 
Fig. 1^ for 3'rd-order vertex and in Fig. |l^ for the 4th order 
vertex) and then to neglect the fully reducible contribu- 
tion (empty objects) for n -\- 1-order vertex in order to 
select that appropriate infinite partial resummation that 
does not clash with the rescaling symmetry (^) that we 
want to preserve. This will be shown and utilized in the 
forthcoming paper jis]] . In particular we will show that 
when we neglect the fully irreducible contributions to the 
4th order vertex we get closed nonlinear equations for the 
triple vertices as shown in Figs|6|3 and 0. By neglecting 
the fully irreducible contributions to the 5th order vertex 
(see Figs. |l^ as and example) we find closed equations 
for the full 3rd and 4th order vertices, and so on. 
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a. What follows from the Wyld equation: 

d ^ }_ d I 2 
-^12 - M7 ^wwv^ 




FIG. 12. Flow chart describing the diagrammatic deriva- 
tion of the first equation in the hierarchy of equations for the 
correlation functions. 



9. APPENDIX A: DISCUSSION OF THE 4TH 
AND THE 5TH ORDER VERTICES 

There are four types of 4th order vertices, as shown 
Fig. 0a. The various contributions to the 4th order 
vertices are grouped according to their topology in three 
classes as shown in Fig. The first class is known as 
"two-eddy reducible" or "skeleton" contributions. Such 
contributions, denoted by joint triangles, can be split into 
two pieces by cutting off one propagator, see panels c 
and d. The second class is known as "irreducible", and 
is denoted as an empty square. The grey square is the 
third class, termed "partly reducible", and made of the 
irreducible vertex and the reducible contributions. The 
need to consider separately the grey vetex stems from its 
appearance in the exact representation of the third order 
vertices, see Fig.^. All the possible types of skeleton con- 
tributions to the 4th order vertices are shown in Figsj^. 
The diagrammatic notation of the skeleton contribution 
is a mnemonic to stress the fact that they consist of two 
dressed triple vertices which were denoted as black tri- 
angles. The joining at the apex hides a bridge that can 
be made from any of the available propagators, either a 
correlator with a wavy line, or a Green's function with 
a straight-wavy line with either orientation. All these 
terms are presented in Fig. 0d. In contrast. Fig. 0c have 
only one type of propagator, with the straight line on 
the right. The two other possibilities have wavy lines 
on the right, and this requires a 3'rd order vertex with 
three wavy tails which is zero identically. Note that the 
line resummed theory does not have bridges made of two 
consecutive propagators with something in between: such 
contributions are "one-eddy reducible" , and all such di- 
agrams have been already resummed in the Dyson line 
resummation. 

As noted, the skeleton diagrams are "two-eddy re- 
ducible" in the sense that they can be split into two pieces 
by cutting across the bridge. Examples of irreducible dia- 
grams for Dn,4:-n which cannot be cut this way are shown 
in Fig. [T^a. There are infinitely many of them. Again, 
this kind of series can be resummed exactly, as is shown 
in Fig. ^3|b. The price is the appearance of an {n + l)'th 
order vertex, 5th order in the present case. This vertex 
is denoted as a grey pentagon. We will call it "partly 
reducible" vertex. Consider now 5th order vertices, for 
example -D4,i, see Fig. |lj. As discussed above we distin- 
guish here the reducible contribution which is explicitly 
shown in Fig .^a and the partly reducible contribution, 
shown as gray pentagon. The need to use a special no- 
tation (grey pentagon) for the partly reducible vertex is 
again motivated by its appearance in the exact represen- 
tatin of the 4th order vertex, Fig.p^. Its further decom- 
position into reducible parts and a fully irreducible term 
(empty pentagon) is shown in Fig .p^. Examples of re- 
ducible (skeleton) contributions to a 5th order vertex are 
shown in Fig. |l^. 

In a forthcoming paper in collaboration with V. Be- 
linicher we will show that the skeleton contributions to 
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a. Series for irreducible part: 




b. Exact representation of Dn,4-n : 




FIG. 13. Panel a: diagrammatic expansion of the quartic 
dressed vertices in terms of triple dressed vertices A, B, and 
C, and propagators. Panel b. Exact representation of the 
irreducible part of DnA-n via the partly reducible 5th order 
vertices shown in Fig 14b. 



a. D4 ^ vertex decomposition : 

2,3 2,3 9^3 

b. (1,2)- irreducible contribution 

2.^3 3 ^2 4,3 5 3 

FIG. 14. Decomposition of the 5th order vertex /)4,i. 

Panel a: decomposition of 1)4,1 into one skeleton contribu- 
tion (first term) and partly reducible 5th order vortex (grey 
pentagon). Panel b: further decomposition of 1)4,1 into re- 
ducible (skeleton) contributions (first three term) and a fully 
irreducible contribution (empty pentagon). 
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Skeleton contributions: 



FIG. 15. Examples of skeleton contributions to the 5tli or- 
der vertex D4 1 . 



the nth order vertices have the same rescaling symmetry 
as the full nth order vertex. This observation means that 
we can achieve a consistent closure by neglecting the ir- 
reducible contributions at any given step of the vertex 
resummation. It is therefore important to observe that 
the skeleton contributions to the n'th order vertex are al- 
ways made of lower order vertices. As shown above, the 
skeleton contributions to the 4th order vertices consist of 
two triple vertices and those for the 5th order vertices 
have contributions of two types. Those with full triple 
vertices connected by two propagator bridges and one full 
triple vertex with one irreducible 4th order vertex con- 
nected by one propagator bridge. We do not present here 
the skeleton contributions to 6'th order vertices. There 
are four types of them: 1) Four black triangles on a line, 
connected via three propagator bridges, 2) Four black 
triangles in a star configuration with one central triangle 
connected via three bridges with three other triangles, 3) 
Two irreducible 4th order vertices (empty squares) con- 
nected by a propagator, and 4) a triple vertex (black 
triangle) connected via a propagator to a 5th order irre- 
ducible vertex. 
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